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Abstract 

A single nucleotide polymorphism (SNP) in the second intron of human TERT (hTERT), rs2736100, acts as a critical factor in 
hTERT synthesis and activation. The rs2736100 SNP was found to be associated with susceptibility to many cancers. 
Recently, inhibition of telomerase and marked telomere shortening were determined to be closely associated with the 
increasing severity of atherosclerosis. The association between the SNP of rs2736100 and the presence of atherosclerosis 
was evaluated in 84 atherosclerosis patients and 257 healthy controls using multivariate logistic regression analyses. The 
proportion of the GG genotype in atherosclerosis patients (17.9%) was significantly higher than in the control group (9.7%). 
Eight variables, including age, gender, cholesterol, high density lipoprotein, homocysteine, total bilirubin, indirect bilirubin, 
and rs2736100 GG genotype, were associated with atherosclerosis with odds ratios of 1.88, 2.1 1, 1.66, 0.23, 1.27, 1.29, 1.53, 
and 1.74, respectively, using multivariate logistic regression analyses. Homozygous GG was demonstrated to be associated 
with the presence of atherosclerosis in our population. 
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Introduction 

A telomere is a region of repetitive nucleotide sequences at the 
end of each chromatid of most eukaryotic organisms that protects 
the end of the chromosome from deterioration or from fusing with 
neighboring chromosomes [1]. Telomere length varies gready 
between species, from approximately hundreds of base pairs in 
yeast to many kilobases in humans, and usually is composed of 
arrays of guanine-rich, six- to eight-base-pair repeats [2,3]. 
Without telomeres, genomes would progressively lose information 
and be truncated after cell division because enzymes that duplicate 
DNA cannot continue their duplication all the way to the end of 
chromosomes [1-3]. Therefore, telomeres act as a disposable 
buffer zone at each end of the chromosomes [4] . Nevertheless, 
telomeres will be consumed during cell division [1-4]. Telomerase 
reverse transcriptase (TERT) is a catalytic subunit of the 
telomerase, which together with the telomerase RNA component 
(TERC), are the most important components of the telomerase 
complex [5]. Telomerase is a ribonucleoprotein polymerase that 
maintains telomere ends by addition of the telomere repeat 
TTAGGG [6] . The enzyme consists of a protein component with 
reverse transcriptase activity, encoded by the gene, and an RNA 
component that serves as a template for the telomere repeat [5,6] . 
Telomerase expression plays a role in cellular senescence because 
it is normally repressed in postnatal somatic cells, resulting in 
progressive shortening of telomeres [5,6]. 



The telomere length of a species is determined by the heredity 
and frequency of cell division [7] . There are theories that claim 
that steady shortening of telomeres with each replication in 
somatic cells may have a role in senescence and in the prevention 
of cancer [7]. Chronic inflammation could accelerate the telomere 
to shorten via increased cell division frequency, and oxidative 
stress could directly result in telomere deletion [8]. 

Human TERT (hTERT) is located in 5p 15.33 [9], and 
rs2736100 is located in the second intron of hTERT. Reports 
have shown that rs2736100 acts as a critical factor in hTERT 
synthesis, in charge of hTERT activation [10]. A single nucleotide 
polymorphism (SNP) of rs2736100 was identified to be associated 
with susceptibility to idiopathic pulmonary fibrosis [11], lung 
cancer [12,13], glioma [14], testicular germ cell cancer, and 
bladder cancer [15,16]. 

Atherosclerosis is the principal cause of myocardial infarction, 
stroke, and peripheral vascular disease, and accounts for nearly 
half of all mortality in developed countries [17—19]. Atheroscle- 
rosis is the buildup of inflammatory cells (macrophages and T 
lymphocytes), vascular smooth muscle cells (VSMCs), and 
intracellular and extracellular lipids on the inner walls of arteries 
[17-19]. These plaques can restrict blood flow to the heart muscle 
by physically clogging the artery or by causing abnormal artery 
tone and function [17-19]. Abnormal VSMC proliferation is 
thought to contribute to the pathogenesis of vascular occlusive 
lesions, including atherosclerosis [20]. Vigorous division and 
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remodeling of VSMCs might consume telomeres. Recently, 
inhibition of telomerase and marked telomere shortening were 
found to be closely associated with the increasing severity of 
atherosclerosis [21]. Because hTERT is the key molecular 
complex that maintains telomere stabilization, genetic polymor- 
phisms in hTERT might be associated with atherosclerosis. In this 
report, the relationship between a SNP of rs2736100 and 
atherosclerosis was evaluated in a case-control study. 

Materials and Methods 

Patients 

To assess genetic polymorphisms related to atherosclerosis, we 
recruited 84 ethnic Han Chinese patients (66 males and 18 
females; age range, 36-79 years; median age, 62 years) with 
atherosclerosis who were unrelated consecutive inpatients at the 
People's Hospital of Yuxi City, China between November 2010 
and November 2012. To obtain an estimate of the genetic 
distribution of the reference allele in the general population in 
China, we also randomly obtained DNA samples from 257 healthy 
individuals who visited People's Hospital of Yuxi City (153 men 
and 104 women; age range, 18-58 years; median age, 34 years) 
with no history of atherosclerosis. The review board of People's 
Hospital of Yuxi City approved this study (approval number: 
YNYXH20 10-003). Written informed consent was obtained 
according to the guidance of the Chinese National Ethics 
Regulation Committee. Participants were simultaneously informed 
of their right to repeal consent by themselves or their kin, 
caretakers, or guardians. 

Patients were diagnosed with atherosclerosis according to 
American Heart Association guidelines [22]. Individuals with 
high blood pressure (Hypertension is defined as systolic pressure 
greater than 1 40 mm Hg or diastolic pressure greater than 90 mm 
Hg.), alcohol abuse, diabetes, and smoking history were rejected 
from the study group. All 84 atherosclerosis patients were further 
confirmed using radionuclide angiography. The 257 healthy 
controls did not have a history of chronic diseases, autoimmune 
diseases, and any cardiovascular diseases. The following clinical 
parameters were obtained for each patient at the time of whole- 
blood collection: age, gender, triglycerides, cholesterol, high 
density lipoprotein, low density lipoprotein, apolipoprotein Al, 
apolipoprotein B, lipoprotein (a), uric acid, homocysteine, total 
bilirubin, y-glutamyltransferase, and indirect bilirubin. 

Polymorphism genotyping 

Genomic DNA was extracted from 100 ju.1 of whole blood using 
the QIAamp DNA Blood Mini Kit (Gaithersburg, MD) according 
to the manufacturer's instructions. Extracted DNA was dissolved 
in 20 ul of 10 mil Tris-HCl buffer (pH 8.0) containing 1 mM 
EDTA, and DNA samples were stored at — 30°C until use. 

Genetic polymorphisms in hTERT were determined by 
MassARRAY®System (San Diego, CA) according to the manu- 
facturer's user guide. We used the primers F (5'-ACGTTG- 
GATGTGACACCCCCACAAGCTAAG-3') and R (5'- 
ACGTTGGATGACAAAGGAGGAAAAGCAGGG-3 ') to am- 
plify the specific hTERT fragment that covers rs2736100 
(Figure 1). A 90-bp hTERT fragment was PCR amplified with 
extracted genomic DNA as the template. PCR was performed 
using HiFifast DNA polymerase (Biovisualab Inc, Shanghai, 
China). The thermocycling conditions were as follows: 94°C for 
10 min, followed by 35 cycles of 94°C for 30 s, 57°C for 60 s, and 
72°C for 60 s. To verify the size of the PCR product, amplicons 
were visualized on 12.5% polyacrylamide gels with appropriate 
size markers. 



The MassARRAY system is based on single-base primer 
extension technology. The MassARRAY technology uses matrix- 
assisted laser desorption ionization time-of-flight (MALDI-TOF) 
mass spectrometry to measure directly the mass of the extension 
product(s) and then correlates the detected mass with a specific 
genotype. The extension primer sequence is 5'- 
TCCGTGTTGAGTGTTTCT-3 ' , and for details of the protocol, 
please refer to SNP Genotyping Using the Sequenom MassAR- 
RAY iPLEX Platform. 

Statistical analyses 

We used a multivariate logistic regression model to calculate the 
statistical power required to detect the contribution of a SNP to 
atherosclerosis risk while including other known risk factors for 
atherosclerosis. SNP status was assigned as X = 0, 1, or 2, which 
represent homozygous for an allele, heterozygous, or homozygous 
for a different allele, respectively. The required sample size, n, for 
the multivariate logistic regression analysis was calculated using 
the following formula [23,24]: 

(z^l+Zj-p) 

[P(l-P)p*2(l-R 2 )] 

where Z u is the upper uth percentile of the standard normal 
distribution; P is the proportion of patients with atherosclerosis 
when X = 1 ; [3* is the size of the effect of an SNP; and R 2 is the 
multiple correlation coefficient (we used 0.009 based on our data) 
relating X with the following covariates: age, male gender, 
cholesterol, high density lipoprotein, homocysteine, total bilirubin 
and indirect bilirubin [25]. When the effect size of a SNP was 
assumed to be 0.69, which corresponded to an odds ratio (OR) of 
2, the required sample size was calculated to be 75 or 100 for a 
statistical power of 80% or 90%, respectively. Based on these 
calculations, our sample size was sufficient for conditions in which 
the OR of an SNP exceeded 2. 

Associations between the clinical parameters (age, gender, 
triglycerides, cholesterol, high density lipoprotein, low density 
lipoprotein, apolipoprotein Al, apolipoprotein B, lipoprotein (a), 
uric acid, homocysteine, total bilirubin, y-glutamyltransferase, 
indirect bilirubin, and rs2736100 genotype) and the presence of 
atherosclerosis were evaluated using Student's t-tests, Mann- 
Whitney U tests, and % 2 tests. Associations between the genotype 
of each locus and the presence of atherosclerosis were evaluated 
using % 2 tests. The Cochran— Armitage test was used to test for 
trends. Possible confounding effects among these variables were 
adjusted using a multivariate logistic regression model, and the 
ORs and 95% confidence intervals (CIs) were calculated. P<0.05 
was considered significant in the two-tailed tests. The Hardy- 
Weinberg equilibria of alleles at individual locus were evaluated 
using the HWE program.3 (ftp://linkage.rockefeller.edu/ 
software). 

Results 

Patients' characteristics 

Of the 84 atherosclerosis patients, 66 were male and 1 8 were 
female, with an age range from 36 to 79 years and a median age of 
62 years. Of the 257 healthy individuals, 153 were male and 104 
were female, with an age range from 18 to 58 years and a median 
age of 34 years. As shown in Table 1, there were no significant 
differences in triglycerides, apolipoprotein Al, or direct bilirubin 
levels between the groups with and without atherosclerosis. In the 
atherosclerosis group, cholesterol, low-density lipoprotein, apoli- 
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Figure 1 . Location of rs27361 00 in the hTERT Gene. Exon and intron organization of the hTERT gene, including the location of rs27361 00. Blue 

rectangles indicate the 16 exons, spaces between two exons were introns. 

doi:10.1371/journal.pone.0085719.g001 



Table 1. Clinical biochemical indexes associated with atherosclerosis. 



Variable 


Total (n = 341) 


Case (n = 84) 


Control (n = 257) 


P 


Age 


40 (18-79) 


62 (36-79) 


34 (18-58) 


0.015 


Gender, male (female) 


219 (122) 


66 (18) 


153 (104) 


0.002 


Triglycerdes (mmol/L) 


2.1 8± 1.60 


1.93 ±1.55 


2.26±1.62 


0.096 


Cholesterol (mmol/L) 


4.41 ±1.06 


4.84 ±1.29 


4.26±0.93 


<0.001 


High density lipoprotein (mmol/L) 


1.08±0.28 


1.17±0.27 


1.04±0.28 


0.001 


Low density lipoprotein (mmol/L) 


2.38±0.73 


2.63 ±0.86 


2.30±0.66 


0.002 


Apolipoprotein A1 (mmol/L) 


1.23 ±0.24 


1.25±0.26 


1.22±0.23 


0.43 


Apolipoprotein B (mmol/L) 


0.69±0.21 


0.83±0.25 


0.64±0.18 


<0.001 


Lipoprotein (a) (mg/L) 


163±164 


246±192 


135±144 


<0.001 


Uric acid (umol/L) 


333±92 


370±85 


321 ±91 


<0.001 


Homocysteine (umol/L) 


15.78±6.36 


21.20±7.68 


14.01 ±4.68 


<0.001 


Total bilirubin (umol/L) 


13.20±5.72 


9.75±3.32 


14.33±5.90 


<0.001 


Direct bilirubin (umol/L) 


4.03 ±1.59 


3.87±1.18 


4.08 ±1.70 


0.219 


y-glutamyltransferase (U/L) 


33.31 ±24.88 


44.69±35.94 


29.6 ±18.63 


<0.001 


Indirect bilirubin (umol/L) 


9.1 7 ±4.45 


5.88±2.48 


10.25±4.42 


<0.001 



Normally distributed data were presented as means ± standard deviation (SD), skewed data were presented as the median (interquartile range). Clinical parameters 
were evaluated to determine an association with the presence of atherosclerosis using the t test or the Mann-Whitney U test. 
doi:1 0.1 371 /journal.pone.008571 9.t001 



poprotein B, lipoprotein (a), uric acid, homocysteine, total 
bilirubin, y-glutamyltransferase, and indirect bilirubin levels were 
higher, whereas high-density lipoprotein was lower in the control 
group, which is consistent with current reports [22,25]. 

rs2736100 allele distribution 

The Hardy-Weinberg equilibria of alleles were evaluated using 
the HWE software (ftp://linkage.rockefeller.edu/software). The 
allele frequency in the atherosclerosis and control groups was in 
Hardy-Weinberg equilibrium (P>0.05), suggesting that allele 
frequencies in the case and control populations remain constant 
in their genetic background. The proportion of the GG, GT, and 



Table 2. The frequence of rs2736100 alleles in case and 
control cohort. 





Genotypes 


Case (n = 84) 


Control (n = 257) 


GG 


15 (17.9) 


25 (9.7) 


GT 


48 (57.1) 


171 (66.5) 


TT 


21 (25.0) 


61 (23.7) 


Proportion of rs2736100 alleles are shown 
doi:1 0.1 371 /journal.pone.008571 9.t002 


as frequency (percentage). 



TT genotypes in atherosclerosis patients were 17.9%, 57.1% and 
25.0%, respectively. The proportion of the GG, GT, and TT 
genotypes in control group were 9.7%, 66.5% and 23.7%, 
respectively. The proportion of the GG genotype in atherosclerosis 
patients (17.9%) was significantly higher than in the control group 
(9.7%) (Table 2). There was no significant difference in the G allele 
frequency between atherosclerosis patients (46.4%) and controls 
(43.0%), suggesting that only homozygous GG was associated with 
atherosclerosis. 

Associations of rs2736100 with atherosclerosis 

The rs2736100 GG homozygote (P= 0.041), cholesterol 
>4.4 mmol/L (P<0.001), low-density lipoprotein >2.38 mmol/ 
L (P= 0.002), apolipoprotein B >0.69 mmol/L (P<0.001), 
lipoprotein (a) >163 mg/L (P<0.001), uric acid >333 umol/L 
(P<0.001), homocysteine ==15.78 umol/L (P<0.001), total biliru- 
bin ==13.20 umol/L (P<0.001), y-glutamyltransferase >33.31 U/ 
L (P<0.001), indirect bilirubin level ==9.17 umol/L (P<0.001), 
and high-density lipoprotein < 1.08 mmol/L (P= 0.001) were 
significandy associated with the presence of atherosclerosis by 
Student's t-tests, Mann-Whitney U tests, and /2. To eliminate any 
possible multicollinearity between or among variables, collinearity 
diagnostics were performed before further statistical analysis, 
cholesterol, low-density lipoprotein and apolipoprotein B dis- 
played collinearity, thus only cholesterol was included in following 
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Table 3. Factors associated with presence of atherosclerosis in multivariate analysis. 



Independent variables 


Odd ratio 


95% CI 


P value 


rs2736100 (GG vs.TT) 


1.74 


1.15-4.63 


0.021 


Age 


1.88 


1.32-4.23 


0.032 


Gender 


2.11 


1.01-4.55 


0.031 


High density lipoprotein (2:1.08 mmol/L vs. <1.08 mmol/L) 


0.23 


0.07-0.66 


0.001 


Homocysteine (2:15.78 umol/L vs. <15.78 umol/L) 


1.27 


1.12-3.31 


0.031 


Total bilirubin (£13.20 umol/L vs. <13.20 umol/L) 


1.29 


1.22-3.43 


0.020 


Indirect bilirubin (29.17 umol/L vs. <9.17 umol/L) 


1.53 


1.20-3.63 


0.039 


Cholesterol (>4.4 mmol/L vs. <4.4 mmol/L) 


1.66 


1.13-4.65 


0.032 



95% CI, 95% confidence interval. 
doi:1 0.1 371 /journal.pone.008571 9.t003 



multivariate logistic regression analyses to avoid multicollmearity 
among cholesterol, low-density lipoprotein and apolipoprotein B. 

To further evaluate the effect of the rs2736100 GG genotype on 
atherosclerosis, multivariate logistic regression analyses were 
performed using these variables (rs2736100, age, geder, choles- 
terol, uric acid, homocysteine, total bilirubin, y-glutamyltransfer- 
ase, indirect bilirubin, and high-density lipoprotein). Eight 
variables, including age, gender, cholesterol, high density lipopro- 
tein, homocysteine, total bilirubin, indirect bilirubin, and 
rs2736100 GG genotype, were included in the final model with 
odds ratios of 1.88 (&50 yrs vs. <50 yrs), 2.1 1 (male vs. <female), 
1.66 (>4.4 mmol/L vs. <4.4 mmol/L), 0.23 (>1.08 mmol/L vs. 
<1.08 mmol/L), 1.27 (>15.78 umol/L vs. <15.78 umol/L), 1.29 
(>13. 20 umol/L vs. < 13.20 umol/L), 1.53 (>9. 17 umol/L vs. 
<9.17 umol/L) and 1.74 (GG vs. TT) (Table 3), respectively. 

Discussion 

Reports have shown that rs2736100 acts as a critical factor in 
hTERT synthesis [1—3]. Normally, telomerase activity in somatic 
cells is silent or low, whereas in tumor cells, particularly lung 
cancer, telomerase is reactivated to extend the telomere and 
hence, tumor cells achieve proliferation [1—6]. Atherosclerosis is 
one of the key causes of coronary heart disease, and VSMC 
proliferation is a typical physiopathological characteristic of 
atherosclerosis [22,25]. Research has shown high telomerase 
activity in neutrophils from unstable coronary plaques [26] . In this 
study, our data show that SNP of rs2736100, together with other 
known factors, were associated with the presence of arterioscle- 
rosis. Our study demonstrated that the GG genotype of rs2736100 
is an independent factor that is associated with the presence of 
arteriosclerosis (OR =1.74 in our model). There were no 
significant differences in allele G frequency between arterioscle- 
rosis patients (46.4%) and controls (43.0%), suggesting that only 
homozygous GG is associated with arteriosclerosis. 

Of the clinical parameters, cholesterol, low-density lipoprotein 
and apolipoprotein B displayed collinearity, this is due to the 
biochemical relation between serum lipids and apolipoproteins. So 
we only input cholesterol into the multivariate logistic regression 
analysis. Cholesterol and high-density lipoprotein were demon- 
strated to be positively and negatively associated with arterioscle- 
rosis which consistent with current reports [27]. Lipoprotein (a), 



uric acid, homocysteine, total bilirubin, y-glutamyltransferase, and 
indirect bilirubin level were significantly higher in the arterioscle- 
rosis group than the control group. However, when these variables 
were input into multivariate logistic regression analysis model, only 
homocysteine, total bilirubin, and indirect bilirubin were included 
in the final model. Lipoprotein (a), y-glutamyltransferase, and uric 
acid were not independent factors in the final model, which may 
be due to the influence of other cofactors. 

SNP assessments are now included in genome-wide association 
studies and replicate models to identify disease susceptibility loci. 
This methodology is currently the best choice, although it is time 
consuming and costly [28]. Investigation of the relationship 
between a certain disease and a SNP screened using molecular 
and cellular physiological background remains a cost-effective 
approach. In this study, we focused on rs2736100 based on many 
previous reports [11-16]. We demonstrated for the first time the 
association between GG genotype of rs2736100 and arterioscle- 
rosis. The relative principle of the molecular function of 
rs2736100, regarding any possible direct and indirect biological 
impact of rs273100 on the presence of arteriosclerosis, needs to be 
investigated in future studies. Despite the limitations of a cross- 
sectional study, our analyses showed a exact effect of the 
rs2736100 GG genotype on the risk of arteriosclerosis. Our 
multivariate model included most of the previously reported 
arteriosclerosis risk factors plus the rs2736100 polymorphism. This 
implies that our results can be generalized to the Han Chinese 
population. The uncertainty of the odds ratios arising from the 
study design might be resolved in a large-scale population based 
study. 
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